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1.  INTRODUCTION 


On  May  1,  1975  a new  program  was  initiated  under  the  sponsorship  of  the 
Office  of  Naval  Research  (Physics  Program)  to  carry  out  a systematic  study 
on  the  piezo-optic  properties  of  liquids  under  high  pressure.  Since  then 
significant  progress  has  been  made  in  this  program  and  in  what  follows,  a 
brief  summary  of  the  work  accomplished  to  date  is  presented. 


2.  WORK  ACCOMPLISHED 

The  technique  of  high  pressure-optical  interferometry  which  has  been 
developed  and  successfully  employed  by  the  author  and  his  students  for 
studies  on  the  piezo-optic  properties  of  solids,  has  now  been  adapted  for 
similar  studies  on  liquids.  Measurements  on  the  variation  of  the  refractive 
Indices  of  a number  of  liquids  at  room  temperature  have  been  completed.  The 
results  obtained  have  already  been  communicated  in  the  form  of  two  short 
articles  to  Physical  Review  Letters  and  another  to  Review  of  Scientific  Instru- 
ments. One  of  the  Physical  Review  Letters  was  also  submitted  to  ONR  as  part 
of  the  Technical  Report  No.  1,  dated  December  22,  1975.  Copies  of  the  other 
two  articles  are  attached  as  appendices  to  this  proposal.  Hence,  in  what 
follows  only  some  of  their  highlights  are  mentioned  and  wherever  the  des- 
cription in  any  of  the  articles  is  too  brief,  additional  data  and  discussion 
thereof  are  presented  here. 

In  brief  it  has  been  possible  to  carry  out  the  optical  interferometric 
measurements  with  a number  of  liquids  such  as  water,  CCIa,  chlorobenzene, 
n-pentane,  n-hexane,  n-heptane,  n-octane,  n-nonane,  and  n-decane,  right  up 
to  their  high  pressure-freezing  points  or  14  kbars  or  whichever  is  lower. 

In  all  cases  the  refractive  indices  of  the  liquids  were  found  to  increase  with 
pressure  with  pronounced  nonlinearity.  However,  this  nonlinearity  is  very 
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much  reduced  when  An,  the  change  in  refractive  index,  is  plotted  as  a function 
of  the  Lagrangian  strain  n.  In  the  case  of  solids  it  was  found  that  the 
change  in  refractive  index  was  exactly  linear  with  Lagrangian  strain  in  the 
entire  range  of  pressures  investigated.  But  in  these  liquids,  at  low 
pressures  when  n is  small.  An  vs.  n is  linear  but  becomes  nonlinear  at  high 
pressure.  Such  a nonlinearity  can  arise  from  one  or  more  of  the  following 
causes:  (i)  unreliable  P-V  data  used  to  evaluate  n;  (ii)  instead  of  Lagrangian 

strain  Hi  the  Eulerian  strain  or  some  other  representation  of  strain  should  be 
used;  (iii)  wrong  equation  of  state  of  liquids  used;  (iv)  intrinsic  nonlinearity. 
In  order  to  narrow  down  these  alternatives,  all  these  piezo-optic  data  on  all 
the  liquids  were  analyzed  by  least  squares  technique  using  all  the  published 
P-V  data  and  the  various  equations  of  state  widely  used  in  the  literature. 

Table  I lists  the  different  Equations  of  State  considered  in  this  work. 

Similarly  for  every  liquid  studied  all  the  published  P-V  data  were  also  con- 
sidered in  the  analysis. 

Firstly,  such  an  analysis  was  helpful  in  eliminating  some  P-V  data  as 
not  being  sufficiently  accurate  or  reliable.  Further,  we  found  in  every  case 
the  range  of  linearity  between  An  and  strain  is  much  larger  when  the  strain  is 
evaluated  in  the  Eulerian  frame  of  reference  than  with  the  Lagrangian  represen- 
tation, regardless  of  the  type  of  equation  of  state  used.  Therefore,  E (=  e) 
appears  to  be  a more  useful  strain  measure  than  n.  Identical  conclusions  on 
the  usefulness  of  E over  n were  deduced  by  Davies^  from  the  data  on  pressure 
dependence  of  the  elastic  moduli  as  well  as  the  Hugoniot  data  on  MgO. 

Tables  II-V  compare  the  degree  of  the  fit  of  An,  to  various  degrees  of 
polynomials  in  E,  for  water,  CCI4,  n-octane  and  chlorobenzene,  respectively, 
for  a number  of  equations  of  state  discussed  in  the  literature.  It  is  seen  that 
a linear  relation  between  An  and  E is  obtained  with  2nd-0rder  Murnaghan  equation 
as  evidenced  by  the  minimum  value  of  the  sum  of  the  squares  of  residuals  as  well 
as  the  standard  error.  The  2nd-0rder  Birch,  Tait  and  Keane  equations  also  give 
essentially  similar  results,  whereas  other  equations  such  as  Bridgman,  lst-Order 
Birch  (also  called  Birch-Murnaghan  equation),  and  lst-Order  Murnaghan  equations 
give  poor  fits  unless  one  uses  higher  order  polynomials.  The  results  on  other 

liquids  are  similar  and  hence  are  not  reported  here. 

2 3 

As  Anderson  , Chhabildas  and  Ruoff  have  shown  that,  of  all  the  equations 
of  state  reported  in  the  literature  only  Keane's  equation  yields  physically  as 
well  as  thermodynamically  meaningful  values  on  extrapolation  to  very  high 


Table  II.  Least-Squares  Fitting  Results  for  Water  Using 
Adams  P-V  Data  to  Compute  E;  An  Measured  at 
25°C  Using  A5893&. 


Sum  of  Squares 
of  Residuals 


Standard 

Error 


Equation 


0.985 


12.75  x 10 


0.979  • 0.100 
0.964  0.674 


4.959 


lst-Order  Murnaghan 


-0.986 

-0.980  -0.088 

-0.991  -0.338 


-3.685 


-0.986 


2nd-0rder  Murnaghan 


-0.979 


-0.094 

0.227 


-0.976 


1.153 


lst-Order  Birch 


-0.985 

-0.979  -0.097 

-0.956  0.952 


18.36  x 10" 
14.98 


7.384 


0.986 


0.096 


0.979 


2nd-0rder  Birch 


0.180 


0.727 


0.977 


-0.986 


1.23  x 10-4 
0.99 
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Table  III. 


Least-Squares  Fitting  Results  for  CC1.  Using 
Mopsik  P-V  Data  to  Compute  E;  An  Measured  at 
25°C  Using  A5893A. 


Equation 

An  = AE  + BE2  + CE3 

Sum  of  Squares 
of  Residuals 

Standard 

Error 

A 

B 

C 

Tait 

- - - - — — 

■ 

■ 

- 1.079  x 102 

11.84  x 10"7 
10.63 
2.41 

1.76  x 10~4 
1.69 
0.82 

: 

lst-Order  Murnaghan 

t:  j 

-1.487 

-1.475 

-1.541 

- 1.276  x 102 

16.00  x 10"7 
14.08 
2.60 

2.05  x 10"4 

1.95 

0.85 

2nd-0rder  Murnaghan 

-1.476 

-1.487 

-1.474 

-0.374 

0.808 

24.217 

3.55  x 10"7 

2.10 

1.69 

• 

lst-Order  Birch 

E 

0.164 

-2.539 

- 55.345 

4.25  x 10“7 
3.97 
1.81 

1.06  x 10'4 

1.04 

0.71 

I 

1.55  x 10‘7 

0.64  x 10"4 
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Table  IV.  Least-Squares  Fitting  Results  for  n-Octane  Using 
Eduljee  P-V  Data  to  Compute  E;  An  Measured  at 
25°C  Using  A5893A. 


Equation 

An 

= AE  + BE2  + CE3 

Sum  of  Squares 

Standard 

B 

c 

of  Residuals 

Error 

Hi 

27.53  x 10'7 

2.53  x 10"4 

Tait 

BSm 

■ 

27.53 

2.56 

EH 

-11.760 

16.00 

1.97 

EH 

98.49  x 10'7 

4.78  x 10"4 

lst-Order  Murnaghan 

BESS 

-0.208 

88.25 

4.58 

-1.415 

-3.186 

-28.354 

21.25 

2.28 

EH 

El 

32.43  x 10"7 

2.75  x 10"4 

2nd-0rder  Murnaghan 

Ell 

32.37 

2.78 

-1.362 

-12.015 

20.34 

2.23 

EH 

91.03  x 10'7 

4.60  x 10-4 

lst-Order  Birch 

0.423 

48.88 

3.41 

Eh 

2.811 

22.738 

5.80 

1.19 

EH 

Hi 

26.21  x 10"7 

2.47  x 10"4 

2nd-0rder  Birch 

EEB 

26.17 

2.49 

-1.358 

-11.250 

15.62 

1.95 

EH 

2.15  x 10'4 

Keane 

RES 

1 

19.75 

2.17 

-1.354 

E 

- 9.835 

11.69 

1.69 

Table  V:  Least-Squares  Fitting  Results  for  Chlorobenzene  Using 

Bridgman  P-V  Data  to  Compute  E;  An  Measured  at  25°C 
Using  X5893A 


Equation 

An  = AE  + BE2  + CE3 

Sum  of  Squares 
of  Residuals 

Standard 

Error 

A 

B 

c 

Tait 

-1.718 

-1.701 

-1.760 

0.321 

-2.508 

-30.101 

52.51  x 10"7 
39.06 
8.37 

3.58  x 10"4 

3.12 

1.46 

lst-Order  Murnaghan 

-1.721 

-1.700 

-1.798 

ED 

-49.700 

11.53  x 10"6 
9.44 
1.08 

5.30  x 10"4 
4.86 
1.66 

2nd-0rder  Murnaghan 

— 

■ 

0.691 

11.98  x 10-7 

10.57 

10.56 

1.71  x 10"4 

1.63 

1.64 

lst-Order  Birch 

-1.714 

-1.700 

-1.697 

0.253 

0.404 

1.608 

10.29  x 10"7 
1.91 
1.82 

1.58  x 10"4 

0.69 

0.68 

2nd-0rder  Birch 

-1.714 

-1.700 

-1.708 

- 4.071 

12.81  x 10‘7 
3.12 
2.56 

1.77  x 10"4 
0.88 
0.81 

Keane 

— 

H 

- 4.372 

13.49  x 10'7 
2.75 
2.11 

1.81  x 10‘4 

0.83 

0.73 

L 
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II 

I I pressures.  Hence,  in  all  further  work  we  propose  to  use  Keane's  equation  of 

state,  instead  of  2nd-0rder  Murnaghan,  2nd-0rder  Birch  and  Tait  equations. 

This  analysis  has  revealed  that  Keane's  equation  gives  excellent  linear  fit 
between  An  and  E for  all  liquids  studied,  irrespective  of  the  nature  of  the 
liquid,  whether  polar  or  nonpolar  and  whether  composed  of  spherically  symmetric 
or  long  chain  molecules.  Thus,  the  present  optical  measurements  support 
similar  conclusions  derived  from  the  ultrasonic  and  Hugoniot  measurements  on 
solids  regarding  the  relative  validity  and  usefulness  of  these  equations  of 
state. 

2.1  . Publications 

o 

. c 

1.  Piezo-optic  Behavior  of  Hater  and  Carbon  Tetrachloride  Under  High 
Pressure:  By  K.  Vedam  and  P.  Limsuwan,  Phys.  Rev.  Letts.  35,  1014 

(1975). 

2.  Piezo-optic  Behavior  and  the  Equation  of  State  of  Liquids:  By  K. 

Vedam  and  P.  Limsuwan,  submitted  for  publication  in  Phys.  Rev.  Letts. 

3.  Optical  Interferometry  in  Liquids  at  High  Pressures  to  14  Kbars:  By 
K.  Vedam  and  P.  Limsuwan,  accepted  for  publication  in  Rev.  Sci. 
Instruments. 

2.2.  Oral  Presentation  at  Scientific  Meetings 

Piezo-optic  Behavior  and  the  Equation  of  State  of  Liquids:  By  P.  Limsuwan 

and  K.  Vedam,  to  be  presented  at  the  ensuing  "Washington  Area  High  Pressure 
Colloquium"  at  the  U.S.  Naval  Academy,  Annapolis,  Maryland,  October  13,  1976. 


APPENDIX  I 

Accepted  for  Publication  in  Review  of  Scientific  Instruments 

Optical  Interferometry  in  Liquids  at  High  Pressures  to  14  Kbars* 

K.  Vedam  and  Pichet  Limsuwan 
Materials  Research  Laboratory  and  Department  of  Physics 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania  16802 

Optical  interferometric  measurements  on  fluids  at  high  pressures  have 

been  reported  in  the  literature  but  the  maximum  pressure  attained  in  these 

measurements  have  been  2.5  kbars  in  the  case  of  gases^  and  only  1.1  kbars 

(21 

in  the  case  of  liquids.'  ' The  authors  have  developed  an  interferometric 

system  for  use  in  liquids  up  to  a maximum  pressure  of  14  kbars,  and  thus 

(31 

have  been  able  to  observe  interesting  nonlinear  piezo-optical  effects.  ' 

A schematic  diagram  of  the  liquid  cell  interferometer  is  shown  in 
Figure  1.  The  interferometer  consists  of  a stainless  steel  container  [see 
No.  1 in  Figure  1]  about  2 cm  diameter  and  1.8  cm  long.  One  end  of  this  con- 
tainer is  connected  to  a polyethylene  bag  [7]  containing  the  liquid  [6]  under 
study,  with  a Buna-N  0-ring  [9]  providing  a tight  seal.  The  other  end  of 
the  container  is  also  sealed  with  another  Buna-N  0-ring  [9]  and  a stainless 
steel  cap  [12].  A vitreous  silica  disc  [13],  cemented  with  Eastman  910 
adhesive^  to  the  stainless  steel  cap  [12]  acts  as  the  optical  window  to  the 
interferometer  and  also  prevents  the  mixing  of  the  experimental  liquid  [6] 
inside  the  interferometer  with  the  fluid  pressure  medium  on  the  outside. 

The  optical  components  of  the  interferometer  consists  of  an  optical 
flat  [5],  a vitreous  silica  plano-convex  lens  [3]  and  a vitreous  silica 
spacer  [4]  between  the  optical  flat  and  the  lens.  A small  hole  [8]  in  the 
base  of  the  stainless  steel  container  [1]  as  well  as  in  the  optical  flat  [5] 
allows  the  space  between  the  lens  and  the  flat  to  be  filled  with  the  liquid 
in  the  polyethylene  bag.  These  components  are  properly  aligned  with  the  help 
of  the  three  stainless  steel  screws  [10]  and  the  three  miniature^  spring 
plungers  [11],  so  that  localized  Newtonian  fringes  can  be  obtained  between 
the  optical  flat  and  the  convex  side  of  the  lens  with  the  spacer  and  the 
liquid  in  between.  To  increase  the  contrast  of  the  fringes  both  the  optical 


flat  and  the  lens  were  partially  silvered  with  a thin  coating  of  aluminum. 

The  spacer  for  this  interferometer  was  fabricated  from  vitreous  silica  but 
before  fabrication,  the  large  block  of  vitreous  silica  itself  was  used  to 
determine  its  second-  and  third-order  elastic  constants.  With  this  informa- 
tion and  a knowledge  of  the  pressure  in  the  system  as  read  by  a precision 
Heise  gauge  or  Manganin  gauge,  the  exact  thickness  of  the  spacer  of  the 
interferometer  at  any  pressure  could  be  evaluated. 

The  entire  interferometer  assembly  could  be  inserted  into  the  14  kbar 
high  pressure  optical  vessel  which  has  been  already  described^  in  connec- 
tion with  similar  measurements  on  the  variation  of  refractive  index  of 
alakli  halides  with  pressure  to  14  kbars.  Figure  2 is  a schematic  diagram 
of  the  high  pressure  optical  vessel  and  the  experimental  arrangement  for 
measuring  the  variation  of  refractive  index  of  liquids  with  pressure. 

Plexol  262  (di-2  ethyl  hexyl  adipate)^  was  used  as  the  pressure  medium. 

As  the  pressure  in  the  system  is  increased,  the  flexibility  of  the  poly- 
ethylene bag  containing  the  fluid  under  study  allows  the  optical  interfero- 
meter to  be  filled  with  the  liquid  at  all  pressures.  The  maximum  amount  of 
liquid  necessary  for  these  measurements  is  only  about  3 c.c. 

A chromel-alumel  thermocouple  coupled  to  a potentiometer  has  also  been 
added  to  the  rear  plug  of  the  high  pressure  vessel.  Since  the  mass  of  the 
optical  pressure  vessel  is  rather  large  and  since  the  rate  of  compression  can 
be  adjusted  and  controlled  to  be  very  slow,  it  has  been  possible  to  maintain 
the  thermal  equilibrium  at  every  stage  to  within  0.1 °C  and  thus  corrections 
for  thermal  fluctuations  is  not  needed.  For  measurements  at  temperatures 
slightly  warmer  than  ambient  temperature,  the  entire  pressure  chamber  can 
be  heated  by  a small  heating  coil  wrapped  around  it. 

As  the  pressure  in  the  system  is  increased,  the  refractive  index  of  the 
liquid  in  the  interferometer  as  well  as  the  thickness  of  the  spacer  are 
altered  with  consequent  shift  of  the  interference  fringes  across  a fiducial 
mark  on  the  optical  flat  of  the  interferometer.  The  change  in  refractive 
index  An  of  the  liquid  could  be  evaluated  from  the  interference  formula 

An  « (pX  - 2n-At)/2t0 


(1) 


where  p is  the  number  of  fringes  shifted.  At  the  change  in  the  thickness  of 
the  vitreous  silica  spacer,  and  t the  initial  thickness  of  the  spacer. 

Such  a computational  procedure  implies  the  following  two  assumptions: 

(i)  the  radius  of  curvature  of  the  lens  does  not  change  significantly  by 
the  change  in  pressure  and  (ii)  the  thickness  of  the  spacer  is  the  same  as 
the  true  thickness  of  the  interferometer  at  the  fiducial  mark.  Since  the 
radius  of  curvature  of  the  plano-convex  lens  of  the  interferometer  is  very 
large  1000  cm)  and  since  the  pressure  change  in  any  one  measurement  and 
the  consequent  computation  involved  using  Equation  (1),  never  exceeded  1000 
psi,  the  above  assumptions  are  considered  to  be  valid  and  the  consequent 
errors  involved  are  considered  to  be  negligible. 

Table  I presents  the  typical  results  obtained  on  two  liquids,  n-hexane 
and  n-pentane.  In  the  case  of  n-hexane  the  measurements  had  to  be  terminated 
at  11.66  kbars  since  it  froze  at  this  pressure  at  25°C.  The  table  gives  the 
values  of  the  refractive  indices  of  the  liquids  at  atmospheric  pressures, 
the  total  number  of  fringes  shifted  and  the  total  observed  change  An  in  the 
refractive  index  of  the  liquid.  Figure  3 shows  graphically  the  observed 
variation  of  the  refractive  indices  of  these  two  liquids  with  pressure  to 
14  kbars. 
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Table  I.  Variation  of  the  Refractive  Index  of  Liquids  with  Pressure 

at  A5461A. 


Temp. 

Max.  Pressure 

Total  No.  of 

Liquids 

(°C) 

n 

(kbar) 

Fringes  Shifted 

An 

n-hexane 

25 

1.3751 

11.66 

822 

0.1566 

n-pentane 

30 

1.3535 

14.00 

941 

0.1809 

Figure  Captions 

Figure  1:  Hicb  Pressure  Liquid  Cell  Interferometer:  [1]  stainless  steel 

container,  [2]  stainless  steel  ring,  [3]  vitreous  silica  plano- 
convex lens,  [4]  vitreous  silica  spacer,  [5]  optical  flat,  [6] 
liquid  under  study,  [7]  polyethylene  bag,  [8]  entrance  hole, 

[9]  0-rings,  [10]  stainless  steel  screw,  [11]  spring  plunger, 
[12]  stainless  steel  cap,  [13]  vitreous  silica  window,  [14] 
teflon  ring,  and  [15]  stainless  steel  plate. 

Figure  2:  Schematic  diagram  of  the  experimental  arrangement  for  measuring 

the  variation  of  refractive  index  of  liquids  with  pressure. 


Figure  3:  Variation  of  the  refractive  indices  of  n-pentane  and  n-hexane 

with  pressure. 
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Piezo-Optic  Behavior  and  the  Equation  of  State  of  Liquids 
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ABSTRACT 

Interferometric  measurements  on  a number  of  liquids  at  high  pressure 
(to  1A  kbars)  show  (i)  first  experimental  proof  that  the  Eulerian  frame- 
indifferent  analogue  strain  E is  a more  useful  strain  measure  than  the 
Lagrangian  strain,  (ii)  Keane's  equation  of  state  best  describes  each  and 
all  liquids  studied,  (iii)  change  in  refractive  index  vs  E is  linear  for  all 
liquids,  even  though  the  volume  strain  involved  is  as  high  as  33*. 

Consider  the  piezo-optic  behavior  of  materials  in  their  transparent 

1 2 

region  of  spectrum.  In  the  case  of  solids  such  as  alkali-halides  , a-quartz  , 

o 

vitreous  silica  etc.,  the  relationship  between  the  change  in  refractive 

index  An  and  pressure  becomes  slightly  nonlinear  at  high  pressures  (say 

above  5 kbars),  but  the  same  data  exhibit  perfect  linear  relationship  between 

An  and  the  Lagrangian  strain  n in  the  entire  range  of  pressures  studied.  In 

the  case  of  liquids  such  as  water  and  CCl^,  as  mentioned  in  a previous 

4 

article  , An  is  grossly  nonlinear  with  pressure,  and  An  vs_ q is  linear  only 
at  strains  less  than  2 or  3%.  Motivated  by  the  linearity  of  An  vs  n for  the 
solids,  one  can  view  the  nonlinearity  of  the  An  vs  n for  the  liquids  as 
possibly  caused  by  (i)  unreliable  P-V  data  used  to  evaluate  r|»  (ii)  use  of  an 
inappropriate  equation  of  state  for  the  liquids  when  extrapolating  and  inter- 
polating literature  P-V  data  to  cover  our  entire  14  kbar  pressure  range. 


(ill)  nonapplicability  of  Lagrangian  strain  as  a strain  measure  at  our  very 
high  strains,  and  (lv)  perhaps  intrinsic  nonlinearity  of  An  with  respect  to 
all  strain  measures  for  liquids.  Investigation  into  these  possibilities 


was  made  with  high  pressure  interferometric  measurements  on  a number  of 
liquids  under  hydrostatic  pressure  and  the  results  and  conclusions  are  pre- 
sented here. 

Thomsen^  and  Davies**’^  have  recently  reviewed  the  significance  and  pro- 
perties of  the  various  representation  of  strains,  in  particular  the 
Lagrangian  (n)  and  the  Eulerian  strain  (e).  In  brief,  for  the  case  of 
liquids  under  hydrostatic  pressure,  we  have 


n " 2 V 


l ; 


n f 


where  Vq  and  V are  the  initial  and  final  specific  volumes.  As  Truesdell  and 


Nolf 


have  shown,  the  Lagrangian  or  the  "material"  strain  is  rotationaily 


invariant,  i.e.,  invariant  under  changes  of  the  frame  of  reference,  whereas 
the  Eulerian  or  "spatial"  strain  does  not  in  general  satisfy  this  invariancy 
criterion.  Hence,  to  overcome  this  limitation  of  e,  Davies**  has  emphasized 
that  the  frame  indifferent  analogue,  E,  of  the  Eulerian  strain  £ should  be 
used  rather  than  £.  However,  for  the  special  case  of  isotropic  bodies  under 
hydrostatic  pressure  E - E and  the  invariancy  condition  is  trivially  satis- 
fied. Thus  for  the  liquids  in  our  hydrostatic  pressure  experiments,  one  can 
conveniently  describe  a frame-invariant  strain  at  any  pressure  in  either  the 
Lagrangian  representation  or  the  Eulerian  frame  - indifferent  analogue 
representation  (which  here  is  identically  equal  to  Eulerian  strain) . These 
strains  can  be  evaluated  by  using  an  appropriate  equation  of  state  to 


extrapolate  and  Interpolate  experimental  P-V  data  to  cover  our  entire 
14  kbar  accessible  pressure  range. 

The  equations  of  state,  for  liquids  and  solids,  that  are  widely  used 
9 10 

In  the  literature  ’ contain  two  or  three  parameters;  bulk  modulus  plus 
first  and  sometimes  second  pressure  derivatives  of  the  bulk  modulus.  Since 
these  constants  cannot  always  be  obtained  with  high  accuracy  from  even  the 

9 

best  available  data  , they  were  evaluated  for  our  studies  by  least  squares 

fitting  each  equation  of  state  to  literature  P-V  data  for  each  liquid.  The 

standard  error  of  estimate  of  each  fitting  was  used  as  an  indication  of  the 

ability  of  each  equation  of  state  to  describe  that  particular  liquid.  In 

addition,  where  more  than  one  source  of  P-V  data  was  available  for  a liquid, 

the  standard  error  was  used  to  rank  the  reliability  of  each  P-V  data  set. 

Table  I presents  some  of  the  experimental  results  obtained  on  a number 

4 11 

of  liquids  by  the  high  pressure  interferometric  method  * . In  every  case 

except  n-pentane  the  maximum  pressures  listed  were  determined  by  the  freezing 

of  the  liquid  at  the  stated  temperature.  It  is  found  that  the  range  of 

linearity  between  An  and  strain  is  much  larger  when  using  the  strain  E instead 

of  ri.  This  is  true  for  each  and  all  equations  of  state  fitted  to  all  P-V 

data  of  all  the  liquids  studied.  (Note  that  we  will  later  show  that  three 

equations  of  state  gave  truly  linear  plots  of  An  xs.  E over  the  entire  14 

kbar  pressure  range  for  all  liquids  studied).  Figure  1 shows  a representative 

12 

graph  of  An  v£  strain,  for  water  evaluated  from  the  P-V  data  of  Adams  and 
the  2nd-0rder  Murnaghan  equation  of  state.  It  is  seen  that  An  increases 
truly  linearly  with  the  strain  parameter  E for  the  entire  stability  field  of 
water  in  the  liquid  phase,  but  that  An  vs.  n is  linear  only  below  about  3% 
strain.  Similar  results  were  obtained  with  every  liquid  studied  thus  far. 


An  important  result  of  this  work  is  that  we  have  experimentally  shown 


i.. 


for  the  first  time  that  E is  a more  useful  strain  measure  than  the 
Lagrangian  strain  r).  Identical  conclusions  for  the  usefulness  of  E over  q 
were  theoretically  deduced  by  Davies^  from  both  the  ultrasonic  data  on  the 
pressure  dependence  of  the  elastic  moduli  and  the  Hugoniot  shock  wave  data 
for  MgO. 

Table  II  compares  the  fit  of  An  for  water  to  various  degrees  of  poly- 
nomials in  E for  six  of  the  widely  used  equations  of  state  discussed  in  the 
literature.  It  is  seen  that  a good  linear  relation  between  An  and  E is 
obtained  with  the  Tait,  2nd-0rder  Murnaghan,  2nd-0rder  Birch,  and  Keane 
equations  of  state  as  evidenced  by  the  value  of  the  sum  of  the  squares  of 
the  residuals  as  well  as  the  standard  error  of  estimate.  The  other  equa- 
tions, including  those  not  listed  there  (such  as  Bridgman  equation,  etc.) 
yield  poor  linear  fit  and  require  higher  order  polynomials  in  E. 

Similar  analysis  with  the  other  liquids  show  that  only  three  equations 

of  state  give  a good  linear  relation  of  An  with  respect  to  E for  all 

liquids.  They  are  the  2nd-0rder  Murnaghan,  2nd-0rder  Birch,  and  the  Keane 

13  14 

equations  of  state.  Anderson  , and  Chhabildas  and  Ruoff  have  shown  that 
of  all  the  equations  of  state  reported  in  the  literature,  only  the  Keane's 
equation  yields  physically  as  well  as  thermodynamically  meaningful  values  on 
extrapolation  to  very  high  pressures.  Hence  the  Keane  equation  is  preferred 
over  the  other  two  that  gave  linear  An  vs_ E plots. 

It  is  one  of  the  more  important  results  of  our  experiments  that  a single 
equation  of  state,  Keane's  equation,  is  found  to  give  an  excellent  linear  fit 
between  An  and  E for  all  the  liquids  studied,  irrespective  of  the  nature  of 
the  liquid,  whether  polar  or  nonpolar  and  whether  composed  of  spherically 
symmetric,  or  planar,  or  long  chain  molecules  and  even  though  the  volume 


I 


strain  involved  is  as  high  as  33%.  Table  II  shows  the  linearity  of  An  vs  E 
calculated  using  the  Keane  equation  of  state. 

Since  one  generally  uses  the  same  set  of  equations  of  state  for  both 
liquids  and  solids  it  is  not  unreasonable  to  expect  the  same  general  con- 
clusions derived  above  for  liquids  to  be  valid  for  solids  under  hydrostatic 
pressure  as  well.  The  results  on  MgO^  and  NaCl^  support  this  conjecture. 
Finally,  the  conclusions  of  the  present  investigation  should  be  particularly 
useful  now  due  to  the  renewed  activities  in  the  physics  of  condensed  matter 

at  high  hydrostatic  pressure  with  the  development  of  diamond  anvil  high 
15 

pressure  system 
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